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ABSTRACT 

We present exact radially self-similar solutions of special-relativistic magnetohydrodynamics representing 
"hot" super-Alfvenic outflows from strongly magnetized, rotating compact objects. We argue that such outflows 
can plausibly arise in gamma-ray burst (GRB) sources and demonstrate that, just as in the case of the trans- 
Alfvenic flows considered in the companion paper, they can attain Lorentz factors that correspond to a rough 
equipartition between the Poynting and kinetic-energy fluxes and become cylindrically collimated on scales com- 
patible with GRB observations. As in the trans-Alfvenic case, the initial acceleration is thermal, but, in contrast 
to the solutions presented in Paper I, part of the enthalpy flux is transformed into Poynting flux during this phase. 
The subsequent, magnetically dominated acceleration can be significantly less rapid than in trans-Alfvenic flows. 

Subject headings: galaxies: jets — gamma rays: bursts — ISM: jets and outflows — MHD — methods: 
analytical — relativity 
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1. INTRODUCTION 

In the companion paper (Vlahakis & Konigl 2003, here- 
after Paper I) we presented a general formulation of special- 
relativistic magnetohydrodynamics (MHD) and derived exact 
radially self-similar solutions for axisymmetric outflows from 
strongly magnetized, rotating compact objects. 1 We pointed 
out that our results should be relevant to relativistic-jet sources 
in a variety of astrophysical settings (including miniquasars and 
active galactic nuclei), but our main focus has been on gamma- 
ray burst (GRB) outflows. 

In Paper I we concentrated on trans-Alfvenic solutions, in 
which the poloidal magnetic field component dominates (or at 
least is not much smaller than) the azimuthal field component 
at the base of the outflow. In the present paper we consider out- 
flows in which \B$/B p \ 1 at the origin. The poloidal com- 
ponent of the Alfven velocity in these flows is very small from 
the start, and consequently they are effectively super-Alfvenic 
throughout. 2 This property of the flows is demonstrated more 
formally in §2. Contopoulos (1995) discussed nonrelativistic, 
"cold" outflows of this type in the case where B p = and 0. 
Here we study super-Alfvenic outflows in the relativistic limit, 
taking into account thermal-pressure effects and allowing for 
the presence of a finite poloidal field component (which is re- 
quired for angular-momentum transport by the flow). 

As in Paper I and in Vlahakis & Konigl (2001), we consider 
a debris disk around a black hole with a massive-star progen- 
itor as a specific example of a long-duration GRB source. If 
the outflow emanates from a region of characteristic cylindrical 
radius 03, and radial width ( ACS), (where the subscript i denotes 
quantities evaluated at the disk surface), and if the disk injects 
energy into the flow mostly in the form of a Poynting flux over 
a time interval Af , then the total injected energy (from the two 
surfaces of the disk) is 



where E = BpV^/c-B^Vp/c is the electric field amplitude. In 
the trans-Alfvenic case E w B p V^ j c, and one infers 
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where the factor (B p i/\B^i\) on the right-hand side is typically 
> 1 for trans-Alfvenic flows. 

A massive star could give rise to a black-hole/debris-disk 
system in at least three ways: (1) a prompt (within ~ Is of 
core collapse) formation in which no supernova shock is pro- 
duced (the "Type I collapsar" scenario; e.g., MacFadyen & 
Woosley 1999); (2) a delayed (within 30- 3000s of core col- 
lapse) formation resulting from fallback in a failed supernova 
explosion (the "Type II collapsar" scenario; e.g., Fryer 1999); 
(3) a strongly delayed (by a period of weeks to years) forma- 
tion resulting from the collapse (following the loss of rotational 
support through a pulsar wind) of a rapidly rotating, massive 
neutron star created in a successful supernova explosion (the 
"supranova" scenario; e.g., Vietri & Stella 1998). In the latter 
scenario, the initial magnitude of B p in the debris disk would 
not exceed typical radio-pulsar values (~ 10 12 - 10 I3 G), and 
in the Type I collapsar picture the magnitude of an inward- 
advected large-scale poloidal field would likely be even lower. 
In the Type II collapsar scenario it is in principle possible for 
a strong poloidal field to develop on the scale of the stalled su- 
pernova shock on account of either convection (e.g., Thompson 
& Murray 2001) or a hydromagnetic instability (e.g., Akiyama 
et al. 2003), but the ability of these mechanisms to produce or- 
dered disk poloidal fields of the required magnitude remains 
questionable. In any case, most of the long-duration GRBs ob- 
served to date are interpreted in the collapsar picture as being 
of the Type I variety (MacFadyen, Woosley, & Heger 2001). 



SittcEiBtitBiiAtyiAt, (1) 

1 In this paper we continue to use the notation of Paper I; we identify equations, sections, and figures from that paper by placing the numeral I in front of their 
respective numbers. 

2 The tendency of the Alfven point in trans-Alfvenic outflows to move close to the disk surface when the value of \B P /B,p\ at the base of the flow decreases below 
~ 1 was illustrated (in the nonrelativistic context) by Cao & Spruit (1994). 
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The most promising mechanism of creating strong fields 
in debris disks is through shear amplification of the seed 
poloidal field by the strong (Keplerian) differential rotation 
(e.g., Meszaros & Rees 1997). This process generates az- 
imuthal fields that can greatly exceed the poloidal field compo- 
nent inside the disk. Kluzniak & Ruderman (1998) and Ruder- 
man, Tao, & Kluzniak (2000), who considered this mechanism 
in the contexts of both a debris disk and a newly formed neu- 
tron star, suggested that the wound-up field would rise buoy- 
antly but not necessarily axisymmetrically and could emerge 
from the disk (or stellar) surface as a flux rope with \B P i\ of or- 
der the amplified \B^,\. A similar scenario has been proposed for 
the formation of active regions in the sun and has been tested by 
numerical simulations (e.g., Brummell, Cline, & Cataneo 2002; 
Abbett & Fisher 2003). Kluzniak & Ruderman (1998) and Ru- 
derman et al. (2000) also suggested that repeated episodes of 
field windup and escape could be the origin of the observed 
intrinsic GRB variability (see § 1.4.1.1). 

Kluzniak & Ruderman (1998) pointed out that an alterna- 
tive outcome of the field winding process might be the escape 
of disconnected (through magnetic reconnection) toroidal flux 
loops. In this case \B^,i/B P i\ could be ^> 1. This is similar to the 
"plasma gun" scenario proposed by Contopoulos (1995). Shi- 
bata, Tajima, & Matsumoto (1990) argued, on the basis of nu- 
merical simulations, that this behavior is likely to occur explo- 
sively and could be manifested by disks in which the pressure 
of the amplified field increases above the thermal pressure (see 
also Haswell, Tajima, & Sakai 1992 and Chakrabarti & D'Silva 
1994). It is interesting to note that the field windup and recon- 
nection cycle underlying these models was also at the basis of 
early models of disk magnetic viscosity (e.g., Eardley & Light- 
man 1975; Coroniti 1981; Stella & Rosner 1984). More recent 
disk-viscosity models have been based on MHD turbulence in- 
duced by the magnetorotational instability (MRI; e.g., Balbus & 
Hawley 1998). The debris disks invoked in GRB source models 
are envisioned to be initially threaded by a subthermal poloidal 
magnetic field, which could make them susceptible to this in- 
stability. One might then wonder whether most of the energy 
liberated in the disk would be viscously dissipated into heat, 
with little left to power Poynting flux-dominated outflows. Al- 
though a definitive answer to this question would entail a fully 
global numerical investigation of the nonlinear development of 
the MRI, existing (radially localized) simulations (e.g., Miller 
& Stone 2000) indicate that disks threaded initially by a weak 
poloidal field with a nonvanishing mean evolve into configura- 
tions that are magnetically dominated throughout (and thus no 
longer MRI-unstable). 

In this paper we explore the possibility that at least some 
GRB outflows are produced with \B^,i/B P i\ 3> 1. In this case 
E w —B^Vp/c, and equation (1) implies 
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where V P i is estimated to be of the order of the speed of sound 
( ~ c/V5) at the base of the flow. Although the behavior of B$ ; 



in the field-windup scenario is inherently nonsteady, it can nev- 
ertheless be incorporated into the quasi-steady "frozen pulse" 
formalism presented in Paper I, as we proved in § 1.2. 1 . 

A brief overview of the application of the r self-similar 
model outlined in Paper I to super-Alfvenic flows is given in 
§ 2. We present representative "hot" and "cold" solutions in § 3 
and discuss analytic approximations to our results in § 4. Our 
conclusions are summarized in § 5. 

2. OUTFLOW DESCRIPTION 

As in Paper I, we consider a multiple-shell outflow con- 
sisting of baryonic matter, electrons (that neutralize the pro- 
tons), e pairs in thermodynamic equilibrium with photons, 
and a large-scale electromagnetic field. We utilize the r self- 
similar, ideal-MHD model described in § 1.3. Adopting as typ- 
ical values w 3 x 10 14 G and B pi w 3 x 10 12 G, we have 
x = E /B p w —B^Vp/ cB p w 60 3> 1 at the base of the flow, so the 
condition M 2 + x 2 > 1 for super-Alfvenic flows (see § 1.3.2.1) 
is satisfied. In fact, for typical parameters the "Alfvenic" Mach 
number M (eq. [1.15]) is also ^> 1 at the base of the flow, which 
further sharpens the above inequality. 

From equation (1.9) = Wfl- V p (—B^/B p ). At the base of a 
trans- Alfvenic flow B p > —B^ and hence O w V^/tS. However, 
when B p <C -B#, the last term on the right-hand side of the ex- 
pression for Vtf, is 3> c, so the first term on the same side must 
also be ^> c in order for their difference to remain < c. 3 This 
implies that w —B ^V p / (B3B p ) ^> V^/OJ at the base of a super- 
Alfvenic flow. In contrast to the trans-Alfvenic situation, in this 
case the fieldline constant £1 cannot be identified with the matter 
angular velocity at the footpoint of the fieldline at the midplane 
of the disk. As was already remarked upon by Contopoulos 
(1995), the fieldline angular velocity il is generally not a mean- 
ingful material property. We prefer instead to interpret it in 
terms of the electric potential $: E = -(Q/c)V s A = -(l/c)V,<I>, 
son = (d/dA)<i>(A,s). 

In super-Alfvenic flows, E w -B^Vp/c w —B$ everywhere 
above the disk. This means that the transfield components of 
the electric and magnetic forces are comparable, f £ j_ w -fgj.. 
This, in turn, is conducive to a successful evolution of the out- 
flow, as the strong magnetic collimation that tends to pinch the 
flow is countered by the decollimating electric force, which al- 
lows the jet to reach large distances. 4 

The expressions for the various physical quantities that ap- 
pear in the r self-similar model can be simplified by taking ac- 
count of the inequalities M 2 3> 1 -x 2 and i>i^, which apply 
in the super-Alfvenic regime. We present the results in Ap- 
pendix A, showing explicitly the dependences on r, 8, and s. 
For example, the expression for B^ has the form 



B, 



,= C 

const 



x r F ~ 2 x 



(sin 6) 



F-2 



X F ~ i {X 2 + M 2 ) 



x[g(ct-£)] 



s=ct-l 



Because of the s dependence, the outflow looks like an outward- 
propagating wave, which is the generalization of the elec- 
tromagnetic outflows (ignoring inertia) described by Lyutikov 
& Blandford (2002). 5 The unknown functions of 9 can be 
obtained after integrating the ordinary differential equations 



3 The implied inequality Q3S1 > c is another manifestation of the fact that the region where B p <C — B^ is located downstream from the "light cylinder," and hence is 
super-Alfvenic. 

4 For Vpi <C c, corresponding to E <g; -B^, we obtain nonrelativistic solutions, which resemble the ones found by Contopoulos (1995). 

5 Their force-free formulation enabled them to also include the oppositely-directed [oc g(ct + £)] wave, which, however, is ■y 2 times weaker than the outgoing one. 
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Table 1 

Parameters of Representative Solutions 1 " 



solution F 


e.O 


e, 


po;(gcnr j ) B pi (G) 


B H (G) 




Vpi/c 




a 1.01 


85 


0.7685 


726 1.2 Xl0 u 


-1.2 X 10 14 


60 


0.7627 


0.5 


b 1.01 


85 


0.011 


437 1.2 xlO 12 


-1.2 x 10 14 


60 


0.7737 


0.5 


T In both cases T 


= 4/3,z c 


= 0, 05, , m 


= 1.5 X 10 b cm, and r ovA /r m = 


3. 
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FIG. 1 . — Main properties of the two representative solutions discussed in §3. See text for details. 
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(I.B2d) and (I.B2e). For a specific solution we need to give 
the value of the model parameter F and seven boundary con- 
ditions at the base of the flow. The latter could be, for ex- 
ample, the values of ©^ponSpM^n^Vpi,^ at the point 
(r, , 6i) where the innermost fieldline is anchored. (The stream- 
line constants p,q,<jM ,Xa, and Bo®l~ F can then be inferred 
from eqs. [I.B7a]-[I.B7c].) The requirement that the outflow 
cross the modified fast-magnetosonic singular surface imposes 
a constraint on these boundary conditions: this is found numer- 
ically as in Paper I. 

3. RESULTS 

We present the results of the numerical integration for two 
representative solutions (labeled a and b), for which the bound- 
ary conditions are given in Table 1 . The most important phys- 
ical quantities are plotted in Figure 1, in which each column 
corresponds to a given solution. The properties of these solu- 
tions are described in detail in the following subsections. 

3.1. Solution a: A Hot, Fast-Rotator Outflow 

This solution is intended to represent a realistic GRB out- 
flow: the total energy initially resides predominantly in the 
electromagnetic field, but the thermal part is nonnegligible, cor- 
responding to temperatures 0, < 1 . 

Figure lal shows the meridional projections of the innermost 
(anchored at G3,- 4„ = 1.5 x 10 6 cm) and outermost (anchored at 
05; .out = 4.5 x 10 6 cm) fieldlines. Over much of the flow the 
line shape is fitted by (z/10 5 cm) w 4(03/ 10 6 cm) 1 ' 48 , and it be- 
comes cylindrical asymptotically. The dashed lines labeled 1 
and 2 represent the optical paths for photons that originate at 
two points on the innermost streamline. The corresponding op- 
tical depths r(l) and r(2) are shown in Figure la5. It is seen 
that the flow becomes optically thin (with r decreasing to ~ 1) 
just before reaching the asymptotic cylindrical regime, corre- 
sponding to maximum baryon loading (see Vlahakis & Konigl 
2001). 

Figure la2 shows the Lorentz factor 7 and the two parts of 
the energy flux in units of 7poc 2 Vp (i.e., the mass-flux xc 2 ): 
the Poynting contribution -EB^ /4irjpocV p and the matter con- 
tribution £7. The ratio of these two parts is the magnetization 
function a = -EB<p j 'Air^ 2 pocV p , whereas their sum is a con- 
stant: fj, = ^—EB^/A-K^pocVp (see eq. [1. 13d]). 6 All quanti- 
ties are given as functions of the cylindrical distance 03 along 
the innermost streamline. It is seen that 7 increases continu- 
ously; for a significant portion of the flow it can be fitted by 
7 i=a (03/10 6 cm) 0,46 , and its asymptotic value is 700 w 500. The 
magnetization function decreases from cr, w 15 near the disk to 
Coo ~ 0.5 asymptotically. 

Figure la3 shows the poloidal and azimuthal components 
of the magnetic field (B p oc 03~ 2 , B^ oc 03" 1 ), the temperature 



(0 oc 03- U Mj ), and the comoving baryon density (po oc 03 ). 

Figure la4 shows the thermal pressure (due to radiation and 
pairs) P, and various electromagnetic pressures. The displayed 
solution corresponds to the case where B 2 /8tt < P < B\ /8tt 
throughout the flow. 

Figure la5 shows the radial profile of the mass-loss rate. For 
03 (jO ut/f<3(,in = '"out/'"in = 3, M k, 2.7 x 10~ 7 M Q s~\ correspond- 
ing to a total ejected baryonic mass of M b w 2.7 x 10" 6 M Q for a 
typical burst duration of Af w 10s, and to a total injected energy 
£i = fiM b c 2 w 3.5 x 10 51 ergs. A fraction (1 + cr^y 1 w 66.7% 

6 This constant may, however be different for different shells, i.e., /i = 



of £i is converted asymptotically into baryon kinetic energy 
Ek = JooMbC 2 w 2.3 x 10 51 ergs of the two oppositely directed 
jets (see eq. [20] below). 



3.1.1. Collimation 
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FIG. 2. — (a) Opening half-angle of the jet as a function of the Lorentz fac- 
tor, (b) The quantity j 2 m/TZ (where TZ is the poloidal curvature radius) as a 
function of G5. Both figures correspond to solution a. 



It is generally much more difficult to achieve collimation in 
the case of relativistic flows than in the nonrelativistic case, 
both because the electric force is as large as the transverse mag- 
netic force and almost cancels it (e.g., Chiueh, Li, & Begelman 
1998; Vlahakis 2003) and because the effective matter inertia 
is larger (e.g., Bogovalov 2001). However, for a flow that is 
ejected with a mildly relativistic velocity, the bulk of the colli- 
mation can take place before it becomes extremely relativistic. 
Figure 2a shows that most of the collimation in the displayed 
solution happens near the disk: the flow starts with ??, = 60° 
and 7,- = 2.4 at 03; ,j n = 1.5 x 10 6 cm, but by 03 w 10 8 cm, where 
7 w 10, it is already collimated to i? 10°. The collimation 
is completed at larger values of 7, and the streamlines eventu- 
ally become cylindrical. Figure 2b shows that for 03 > 10 8 cm 
(corresponding to 7 > 10) and up to 03 « 10 n cm (where the 
cylindrical regime commences), the quantity 7 2 05/7£. is approx- 
imately constant (Chiueh et al. 1998; see also § 4). Despite the 
large value of the curvature radius, it is still possible for the 
opening half-angle to vanish asymptotically. 

The collimation is the result of the interplay between the 
magnetic and electric forces in the transfield direction. Fig- 
ure 3b shows all the force components in that direction. It is 
seen that both t E ± and t B ± change sign (at 03 w 5 x 10 6 cm and 
7 x 10 6 cm, respectively). In a small region very close the origin 
{e± > (corresponding to a positive charge density, J° > 0) and 
ffl_i_ < (corresponding to the return-current regime /|| > 0). In 
this region the electric force collimates and the magnetic force 
acts to decollimate the flow. At larger distances f £ j_ < (the 
charge density becomes negative, J° < 0) and fg± > (cor- 
responding to the current-carrying regime < 0). Beyond 
03 Rj 10 8 cm these two forces become almost equal (their dif- 
ference is the much smaller force — f/j_), and the collimation 
continues slowly. 
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FIG. 3. — Force densities (multiplied by rj5 3 ) along the poloidal flow (top 
curve) and in the transfield direction (bottom curve), as functions of 03 along 
the innermost streamline of solution a. The various force components are de- 
fined in eq. (1.18). 

3.1.2. Acceleration 

The poloidal force components along the flow are depicted in 
Figure 3 a. Referring also to Figure la2, we identify three dif- 
ferent regimes (as in the corresponding trans-Alfvenic solution 
described in § 1.4.1): 

1 . Thermal acceleration region 

For 03 < 3 x 10 8 cm, the thermal force exceeds the magnetic 
force, f r || > f B ||. The acceleration is thermal: the increase in 7 
is caused by the decrease in the specific enthalpy £. 7 
A new feature of this solution (not encountered in the corre- 
sponding trans-Alfvenic solution presented in Paper I) is that 
part of the enthalpy flux in this regime is transformed into 
Poynting flux: for 05 < 7 x 10 7 cm the { B n component is neg- 
ative (corresponding to J± < 0), resulting in an increase in the 
Poynting flux and a decrease in the enthalpy flux (see Fig. Ia2). 
(In the corresponding regime of solution a of Paper I, the en- 
thalpy flux oc £7 remains constant.) The magnetization function 
increases from a w 15 near the disk to a w 50 at 03 w 7 x 10 7 cm, 
before starting to decrease in the regime where f B M > 0. 
Combining the conservation relations for the total specific an- 
gular momentum and the total specific energy (eqs. [1. 13c] and 
[I.13d]),weget 



^05V^ 

c £ 



(4) 



Note that for the adopted fiducial values the quantity 

xM 2 (Y±_?L 
M 2 +x 2 -xV 4> /c \ c E 



xl-l 



initial opening half-angle i?,-, so that the region where f B || < 
would be more extended in this case. 

2. Magnetic acceleration region 

From the end of the thermal acceleration zone, where £ w 1, up 
to 03 w 4 x 10 11 cm, it is seen from Figure la2 that the Lorentz 
factor continues to increase. The acceleration in this regime is 
due to magnetic effects: Poynting flux is transformed into ki- 
netic energy flux. Figure 3a shows that the force -fc\\ (which 
measures the increase in 7) is equal to the magnetic force fg|| 
(which derives from the decrease in | 05Z?^ |). 
Even though 03/7£. < I/7 2 (see Fig. 2b), in agreement with the 
analysis of Chiueh et al. (1998), their argument that no signifi- 
cant acceleration can take place at large distances from the light 
cylinder is evidently circumvented in this case: the a function 
in the magnetic acceleration region decreases from a w 50 to 
the asymptotic value k> 0.5. This behavior is explained in 
Vlahakis (2003), who analyzes the efficiency of magnetic ac- 
celeration in relation to the fieldline shape and the role of the 
centrifugal force in the transfield force-balance equation. As 
the flow becomes progressively more matter-dominated, an im- 
portant transition is reached (at a = <r c ) when the centrifugal 
force becomes equal to the electromagnetic force component 
in the transfield direction. The acceleration could in principle 
continue beyond that point only if a transition from a positive to 
a negative poloidal curvature were possible. This cannot hap- 
pen in r self-similar models, which therefore have < a c . 
Large asymptotic cylindrical radii correspond to a compara- 
tively smaller centrifugal force and therefore, (on account of 
ffi± + ffi_L = — fc_i_) to a less magnetized asymptotic flow. This 
explains why the super- Alfvenic solutions presented in this pa- 
per attain smaller values of than the corresponding solu- 
tions of Paper I, which generally collimate more efficiently than 
super- Alfvenic flows. 

3. Asymptotic cylindrical region 

At the end of the magnetic acceleration region the flow becomes 
cylindrical: Figure lal shows that the fieldlines converge to a 
constant value of 03 and that subsequently all the flow quantities 
remain constant. 

3.1.3. Time-Dependent Effects 

The disk activity that determines the time variability of a 
GRB remains an open question, so one is not yet in a position to 
accurately specify the time profile of the boundary conditions 
at the base of the flow. We nevertheless present in this subsec- 
tion an illustrative example of how one may recover the time 
dependence imprinted on the flow by the s dependence of the 
boundary conditions at the origin. In particular, we show how 
to obtain the function 7^ , s) [or, using s = ct-£, the function 
7(^,0; see also §1.4.1.1]. 

Equations (Aid) and (Alf) give 



7 (*,s) = 7 (Oo) 



and 



= gO) qjsp) £(l,sp) 
[&e,s)-lf ~ g(s ) q(s) [ftf,s )-l] 3 



(5) 



(6) 



is positive (in contrast to the trans-Alfvenic case, where x A < 1). 
Thus, equation (4) implies that £7 decreases with increasing 
05V^. This explains the behavior of £7 near the base of the flow 
and indicates when this effect would be most pronounced: it 
may be expected that OSV^ would increase faster for a larger 

7 There is a very small region near the origin where -fg is small, meaning that the Lorentz factor remains constant (see also Fig. Ia2). However, even though V 
remains constant, the poloidal velocity increases (while decreases). 



For given g(s), g(s)//j,(s), q(s) and a known solution 
so) so)} for the reference shell so, we can solve this 
system of equations to obtain {£(£, s), 7(^,5')} for all the other 
shells of the outflow. 
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As a concrete example, we choose the functions g(s) and 
g(s) I 'fj,(s) as shown in the inset of Figure 4. In the context of 
the internal-shock model of GRBs, these functions can be in- 
terpreted as the envelope profiles of a series of discrete pulses. 
Since the Poynting flux is oc g(s), we choose g max /^min = 10 4 , 
corresponding to an azimuthal magnetic field amplification fac- 
tor B$ max /B0 min = 10 2 . The profile of g(s)//j,(s) oc jp Q is sim- 
ilar to the one adopted in Paper I. We also set q(s) oc g(s)/j, 2 (s), 
for which choice the solution of equations (5) and (6) near the 
origin (where £ » 1) is j(£ w 0,s) w -f(£ w 0,so). 8 The result- 
ing Lorentz factor profile j(£ , t) is shown in Figure 4. 



1000 



100 




The asymptotic baryon kinetic energy is the same as in so- 
lution a, and the baryon loading is similarly the maximum 
allowed under the condition that the flow be optically thin 
in the asymptotic regime (tqo w otGJooPo oo/m p w 1). Since 

Ek oc po ooT^^L^f ^007^^00, the smaller asymptotic jet 
radius implies a larger value of (and hence a smaller bary- 
onic mass). The actual numbers for solution b are M w 1.75 x 
10~ 7 M Q s" 1 , w 730 (corresponding to (T^ w 0.58), and an 
acceleration efficiency ^/^^(l + iToo)" 1 f«63.3%. In this case 
P < Bp/87r < B0/8tt throughout the flow. 

Figure lb2 shows a slight initial decrease in £7 w 7. This 
behavior is similar to the one exhibited by solution a and again 
corresponds to a negative value of f B ii . In this case, however, the 
Poynting flux increases not at the expense of the thermal energy 
but rather because a fraction of the kinetic energy of rotational 
motion is transformed into electromagnetic energy (with the re- 
mainder of the rotational energy accounting for the increase in 
the kinetic energy of poloidal motion in this region). 

4. DISCUSSION 

The outflow solutions depicted in Figure 1 exhibit a power- 
law behavior over most of their extent. We now show that these 
power-law scalings can be obtained analytically. 

The Bernoulli equation implies V p w c, or, using equation 
(Ale), 

Fa„ sing X 2 + M\ 
jU sin«?-i?) X 2 ~ ' ' 

Well above the disk surface (where 6 <C tt/2 and £ w 1), 
equation (7) yields 



/rflnCjN _ 1? _ , Fa M / 1 
\ c/lnz ) A fi \ a 



(8) 



where we approximated the magnetization function a = 
-(cEB0/47r)/(£7 2 p o c 2 Vp) by A 2 /.M 2 . In the Poynting flux- 
dominated regime a > 1, so the fieldline shape is 



z oc 03 v- f °m . (9) 

The opening half-angle $ oc cj- F,T «AM-f ^m) ) an( j jjjg curV ature 
radius satisfies 
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03 <9 2 C5(z,A) ( B, 



Fa M ( Fa M 



-T- (10) 



The transfield force-balance equation can be written as 



1 io 4 10 5 10 6 10 7 10 8 10 9 10 10 io 11 10 12 10 13 10 14 10 15 

l(cm) 

FIG. 4. — Lorentz factor plotted as a function of the arclength along the 
innermost fieldline of solution a for several values of time since the start of 
the burst. The dashed line represents the "time independent" solution for the 
reference shell sq = 0.9 X 10" cm. The assumed forms of the functions g(s) 
and g(.v) / fi(s) across the ejected "pancake" are shown in the inset. 

For t = 0.5 s, only the shells with s <E [0, ct = 1.5 x 10 10 cm] 
have been ejected. The ejecta occupy the region from I = 
to £ = ct = 1.5 x 10 10 cm, and the first shells have already been 
accelerated. 

At time f = lis, the whole "pancake" s e 
[0, cAf = 3 x 10 n cm] can be seen (assuming a typi- 
cal duration of At = 10s). It occupies the region £ E 
[ct-cAt = 3 x 10 10 cm,cf = 3.3 x 10 n cm]. 

The acceleration continues and at t = 10 4 s the asymptotic 
stage is reached. Although the shells were ejected with the 
same 7,, they now span a wide range of Lorentz factors [cor- 
responding to the linear dependence 700 (^) = ((Too + 1) ' 
which would have the effect of enhancing the emission effi- 
ciency of the resulting internal shocks (e.g., Kobayashi, Piran, 
& Sari 1997). 

3.2. Solution b: The Cold Case w 1) 

The cold solution presented in the second column of Figure 
1 is more magnetized than solution a: the initial value of its 
magnetization function is 07 w 450. As a result, it collimates 
faster and the asymptotic jet radius C3oo is smaller than the cor- 
responding radius in the fiducial solution. 

8 For these forms of the s dependence, we verified that the neglected terms in the momentum equation (I.12e) are very small compared with the other terms of that 
equation. 



C5V4 m 4 

iw v * ln| — 
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7 2j%(M 2 + x 2 ) 
TIB 2 



(11) 



(e.g., Chiueh et al. 1991; Okamoto 2002; Vlahakis 2003). 
Using x 2 /M 2 rts a » 1, -B^ w E = xB p , and assuming that 
-®B<f, /7 = A m T{£) (a form verified by numerical integration at 
least in the current-carrying regime), we get 

05 1 VA I 



2 03 
7 — =m- 



(12) 



From the definition of the fieldline constant a M (= AVl 2 / ^ A c z ) 
and the fact that the flow is Poynting flux-dominated (p, w 
-XSQ.B^/^ac 2 ) we obtain 

A a M 
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so equation (12) yields 



2 03 fi 
7 — = m — 
1Z om 



(14) 



Figure 2b verifies that the right-hand side is a constant (« 1). 
Equation (10) then implies 



7 : 



Fo M 



mF 
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i.e., 



7 oc 03 
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Fo M 



and 7 cx z 



.(Fctm/m) 



(15) 
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fj,-Fa M 

The scaling 7 cx z/05 can also be derived from kinematic con- 
siderations: c/03 - V m dt ~ (c 2 -K 2 ) 1/2 c/f , and for V ra < V z we get 
7 ~ dzjdXSS cx z/05 for a fieldline shape of the form of equation 
(9). 

The solutions derived in §3 have (3 w 0.46 and 0.49 in the 
"hot" and "cold" cases, respectively. The other scalings are 
easier to derive: using equation (13) we get 

/jA£1 



Z^05" 2 , 



-xB„ 



CO M 



-05" 



(17) 



Note that it is, in fact, the deviation from the last two scaling 
relations that gives rise to the acceleration: 

uAQ, AS1 a 

— — + 05^0 = 7cx03 /5 (18) 

CO M CO M 

(where we used the definition [1. 13d] of /j, in the limit £ ~ 1). 
Finally, from the constancy of the mass-to-magnetic flux ratio 
and from the polytropic relation we infer 



7 A) - 



-03- 2 , 9cxp ( ; /3 cx03- ( ' 3+2) / 3 



(19) 



Alt CO M 

The fiducial "hot" trans-Alfvenic solution presented in Paper I 
also exhibits a power-law behavior. Its power-law exponents 
can be obtained as a special case (f3 = 1) of the preceding re- 
sults, corresponding to a parabolic fieldline shape (z cx OS 2 ) and 
7 cx 03 (see Vlahakis & Konigl 2001). As demonstrated above, 
our fiducial super-Alfvenic flow accelerates significantly more 
slowly with distance from the source. 

Equation (8) shows that, as the flow accelerates and the value 
of the magnetization function declines, the jet opening half- 
angle decreases faster than 9. The cylindrical regime is thus 
reached at a finite distance from the source, corresponding to 

Fo M /j, 

7oo = z = H~Fo M , 



Coo = 



l + cr n 



(20) 



H-Fom 

which yield an acceleration efficiency = 7oo/m = (1 + 

Coo)" 1 = l—Fau/fJ,. Higher efficiencies thus correspond to 
larger values of /i/Fom- We can use equation (Ale) to re- 
late this parameter combination to the boundary conditions at 
the base of the flow: h/Fom = csin0,-/^,,-sin(0i — #,•), or, for 
Oi w n/2, [ijFo M ~ c/V-,j. The poloidal velocity cannot be 
very small (V p /c w -E/B^ should be close to 1 to allow the 
flow to reach large distances; see § 2), so a high efficiency re- 
quires a large initial opening half-angle However, a large 
value of ii/Fom also implies a small value of the exponent (3 = 
dln^f j din®, and since 700 is always close to (1, 05^ w 03,/i 1 /^ 
would have to be large. Since z cx 03 1+/3 , this places a practical 
upper limit on it cannot be too large or else the accelera- 
tion would be completed on scales much larger than typically 
inferred for GRB outflows (~ 10 14 - 10 15 cm). 



Drenkhahn & Spruit (2002) examined a pure-B^, strictly ra- 
dial flow without considering the transfield force-balance equa- 
tion. They showed that magnetic energy dissipation (mod- 
eled in a parameterized manner) results in efficient accelera- 
tion 7 cx r° 25 and a final Poynting-to-kinetic energy conversion 
efficiency w 54%. Our fiducial solution exhibits a faster accel- 
eration 7 cx r° 31 and a higher efficiency w 67%. These results 
indicate that both dissipation and the fieldline shape may play 
an important role in the flow acceleration. 



5. CONCLUSION 

Using the radially self-similar relativistic MHD model pre- 
sented in Paper I, we constructed representative "hot" and 
"cold" fast-rotator solutions of super-Alfvenic flows. We ar- 
gued that super-Alfvenic outflows in which \B^,/B P \ >• 1 al- 
ready at the source could plausibly arise in GRB source models 
that invoke differential rotation in a disk or a star to account for 
the large magnetic field amplitudes that are required for consis- 
tency with the observations. We demonstrated that our "frozen 
pulse" formulation, despite being quasi steady, can still capture 
the expected time-variability of the azimuthal field component 
at the source. We showed that for typical source parameters 
the flows convert Poynting flux to kinetic-energy flux with high 
efficiency (w 67% for our fiducial solution) and collimate to 
cylindrical structures on scales > 10 14 cm. Our solutions con- 
firm (some previous statements in the literature notwithstand- 
ing) that significant magnetic acceleration and collimation can 
take place at large distances from the source (well beyond the 
light cylinder), although most of the collimation is achieved be- 
fore the flow becomes extremely relativistic. 

The super-Alfvenic solutions derived in this paper are distin- 
guished from the trans-Alfvenic ones obtained in paper I in two 
main respects: (1) During the initial thermal-acceleration phase 
of the "hot" solution, some of the internal energy is transformed 
into electromagnetic energy even as another part is used to in- 
crease the poloidal velocity. (The same behavior is exhibited by 
the "cold" solution, except that in that case the energy reservoir 
is the bulk initial rotation rather than the thermal energy.) (2) 
During the subsequent magnetic-acceleration phase, the rate of 
increase of the Lorentz factor with distance from the source can 
be significantly lower than in the trans-Alfvenic case; the rate 
of increase of the jet radius with distance is correspondingly 
higher. Overall, however,, the two types of solution are quite 
similar, and we derived analytic scaling relations that describe 
them both (see also Vlahakis & Konigl 2001). We conclude that 
source configurations with either \B p /B^,\ > 1 or \B P /B^\ <C 1 
could in principle produce viable GRB outflows. Another po- 
tentially important aspect of such outflows, namely, the possi- 
bility that their initial composition is highly neutron-rich, can 
also be modeled within the theoretical framework that we have 
constructed. This is discussed in a separate publication (Vla- 
hakis at al. 2003). 
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APPENDIX 
A. PHYSICAL QUANTITIES 

In the super- Alfvenic regime (M 2 » 1 -x 2 , x » x A ) of the r self-similar model described in Paper I, the expressions for the various 
physical quantities as functions of r, 9, and s = ct-£ take the form: 



B = 



g F / 2 (s Q ) 
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X F rsm(9-i)) 



(z cos $ + rosin i9) - 
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ctm(so) om (j) ^c 2 ^ 2 A" 2 ^- 2 )( A" 2 + M 2 ) ' 



(Ala) 
(Alb) 
(Ale) 
(Aid) 
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where x 2 = ;f 2 g(.v), M 2 = M 2 s(,), ^g^' = J °^y M , = ^W^, and the functions ^(j), and <?(,) 

can be specified freely. The functions X 2 (8), M 2 (9), and z?(6*) can be found from the integration of the ordinary differential equations 
(I.B2d) and (I.B2e) of Paper I, whereas £ is given from 



(£-D T 



q(s)~ 



-M 



(Alf) 
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